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ABSTRACT 

We compute the spectral distortions of the Cosmic Microwave Background (CMB) 
arising during the epoch of cosmological hydrogen recombination within the standard 
cosmological (concordance) model for frequencies in the range 1 GHz — 3500 GHz. We 
follow the evolution of the populations of the hydrogen levels including states up to 
principle quantum number n = 30 in the redshift range 500 $J z $J 3500. All angular 
momentum sub-states are treated individually, resulting in a total number of 465 hy- 
drogen levels. The evolution of the matter temperature and the fraction of electrons 
coming from Hell are also included. We present a detailed discussion of the distortions 
arising from the main dipolar transitions, e.g. Lyman and Balmer series, as well as 
the emission due to the two-photon decay of the hydrogen 2s level. Furthermore, we 
investigate the robusteness of the results against changes in the number of shells con- 
sidered. The resulting spectral distortions have a characteristic oscillatory behaviour, 
which might allow experimentalists to separate them from other backgrounds. The 
relative distortion of the spectrum exceeds a value of 10 -7 at wavelengths longer 
than 21 cm. Our results also show the importance of detailed follow-up of the angular 
momentum sub-states, and their effect on the amplitude of the lines. The effect on 
the residual electron fraction is only moderate, and mainly occurs at low redshifts. 
The CMB angular power spectrum is changed by less than 1%. Finally, our compu- 
tations show that if the primordial radiation held is described by a pure blackbody, 
then there is no significant emission from any hydrogen transition at redshifts greater 
than z ~ 2000. This is in contrast to some earlier works, where the existence of a 
'pre-recombination' peak was claimed. 

Key words: cosmic microwave background - cosmology: early Universe - cosmology: 
theory - atomic procceses 



1 INTRODUCTION 

In the last twenty years, many experiments have been 
devoted to the study of the angular fluctuations of the 
Cosmic Microwave Background (CMB) temperature and 
polar ization. Nowadays, the r esults of the WMAP satel- 
lite iSpergel et alJ I2003L I2006T) in conjunction with other 
datasets (observations of type la supernovae or galaxy sur- 
veys) give us the most precise determinations of the cos- 
mological parameters describing our Universe, and consti- 
tute one of the most important confirmations of the hot Big 
Bang theory. One of the predictions of this theory is the ex- 
istence of small spectral distortions of the CMB blackbody 
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spectrum arising as a consequence of the non-equilibrium 
conditions occurring during the epoch of th e primordial hy- 
drogen recombinati on at redshift z ~ 1100 ( Zcld ovich et alJ 
1968: lPeebles! IT968). These authors showed that, because in 
the Wien part of the CMB spectrum the brightness is ex- 
tremely low, the strongest distortions are connected with the 
Lya line (which today should be found at A ~ 170 fjxo), and 
with the 2s level two-photon decay emission (which should 
peak around A ~ 200 fim). There should also be similar 
signatures arising from the recombination of helium (which 
occurs around z ~ 6000 for Helll and z ~ 2500 for Hell), 
but these appear at longer wavelengths than the hydro- 
gen L ya line (e.g. lDubrovich fc Stolvarovlll997l : IWong et alJ 
2006) and they lie one order of magnitude below the hydro- 
gen spectrum, just due to the lower abundance of helium. 



2 J. A. Rubino- Martin, J. Chluba & R.A. Sunyaev 



iDubrovichl lll975t) proposed to look for the spectral dis- 
tortions due to transitions between higher levels of the hy- 
drogen atom (Balmer and higher series), and after this, 
many papers were devoted to the problem of the formation 
of the hydrogen recombinat i on lines in the early Universe 
jLiuha.rskii fc Snnva.evllT9Sa IpThicki fc dell'Anto'nkl llQfll 



(Liubarskii & SunjagvJ WMLRvbicki & doll' Antonio! .1993: 
Dubrovich & Stolvarov 1995, 1997; Boschan & Biltzinger 



19981: iBurginl 120031: IDubrovich fc Shakhvorostoval 12004 
Kholupenko et al.ll2005HWong et al.ll2006l) . Although there 



were contradictory results in the amplitudes and shapes of 
the features, these computations showed that the ampli- 
tude of the spectral distortions arising from transitions of 
the higher series is very small, with typical values around 
~ 10 -6 — 10 -8 . The important point for us is that these dis- 
tortions should have survived until today and might become 
observable in the near future. 

The measurements from the F ar-InfraRed Abso- 
lute Spectrophotomete r (FIRAS, e.g. iFixsen et alJ Il996t 
IFixsen fc Mathedlioolj) on board the COBE satellite indi- 
cate that the CMB frequency spectrum is described by a 
To = 2.725 ± 0.001 K blackbody spectrum with high pre- 
cision. The limits on the Bose-Einstein and Compton dis- 
tortions are \n\ < 9 x 10" 5 (95% CL) and \y\ < 15 x 10~ 6 
(95% CL), respectively. Moreover. IFixsen fcMatheil (120021) 
point out that recent technological progress permits one to 
improve the limits by nearly 2 orders of magnitude. There- 
fore, we are approaching the sensitivities required for detect- 
ing global features in the spectrum, so detecting individual 
ones could be possible in the next future. Therefore, pre- 
cise computations of the shapes, positions and amplitudes 
of these lines are needed. We also mention that there have 
been attempts to measure the CMB spectrum at centimeter 
wave l engths using the A RCADE experiment jKogut et alJ 
120041: IFixsen et aljliooll . 

In most of the existing calculations simplifying as- 
sumptions have been used in order to make the prob- 
lem computationally tractable. For the study of the Lya 



three-level atom is usually adopted (e.g. 


Zcldovich et al. 


1968 


Peebles 1968; Boschan fc Biltzinger 


1998; Wong et al. 


2006) 


. For transitions in the Balmer 


and 


higher se- 



ries, two assumptions are usually ado 


ated (e.g. iBurgin 


2003 


; Dubrovich & Shakhvorostova 2004; 


KholuDenko et al 


2005 


). First, a quasi-static evolution of the populations for 



all levels above n = 2 is assumed, which reduces the problem 
of solving a system of stiff coupled ordinary differential equa- 
tions to an algebraic system. And second, given that one of 
the main difficulties is treating all the energetically degener- 
ate angular momentum sublevels within each shell, only the 
total population of a given shell is computed and statistical 
equilibrium 1 (SE) among the sublevels is assumed. To date, 
th e only attempt to avoid these t wo assumptions was made 
bv lRvbicki fc dell'Antoniol (Il993l) . who considered the case 
of 10-shells. However, their results were not conclusive and 
were affected by numerical uncertainties. In particular until 
now the question whether the assumption of SE is justified 
has not been adressed satisfactory. 



1 This assumption implies that the population of a level (n, I) is 
given by N n i = (21 + l)N n /n 2 , where N n is the total population 
of the shell with principle quantum number n 



It is important to mention that all the codes which com- 
pute the power spectrum of the angular fluctuations of the 
CMB make use of the r ecombination histo ry as obtained 
with the Recfast code iSeager et alJll999i) . These calcu- 
lations are based on a 300-level hydro gen atom in whic h 
SE among the I sublevels was assumed dSeager et aljEoobT) . 
Therefore, it would also be important to understand the 
validity of this approximation and its effects on the recom- 
bination history. 

In our paper, we have obtained the spectral distortions 
from all lines of hydrogen including 30-shells, following the 
evolution of all angular momentum sub-states separately. 
Since we are until now not including collis ions, from an 
estimate following IPengellv fc Seatonl lll964h this seems a 
natural choice in order to keep the calculations simple. Our 
Universe has extremely high entropy, i.e. photons outnum- 
ber the baryons by a factor of ~ 1.6 x 10 9 . This is why 
collisions in general are not important. However, those colli- 
sions connecting different, energetically degenerate angular 
momentum substates within a given shell should play a role 
for sufficiently high shells, but the results presented here 
should remain unaffected. A detailed discussion of this sub- 
ject will be done in a subsequent paper (Chluba et al. 2006, 
in preparation). 

Our computations have permitted us to show for the 
first time that the deviations in the populations of the an- 
gular momentum substates within a given shell from SE are 
important when studying the spectral distortions to the cos- 
mological blackbody, and that the impact on the residual 
ionization fraction of the hydrogen atoms after recombina- 
tion is of the order of a few percent. However, this change is 
sufficiently important to produce modifications of the order 
of few percent on the predicted angular power spectrum, so 
this effect should be taken into account for future experi- 
ments like PLANCK. 

The main results of this paper are summarised in Fig. Q 
and Fig. [5] The hydrogen lines from the epoch of cosmologi- 
cal recombination appear as broad (Ais/v ~ 0.30 — 0.40) fea- 
tures, having high contrast for the Lyman, Balmer, Paschen, 
Brackett and Pfund series. Altogether they produce a spec- 
tral distortion with a characteristic oscillatory behavior, 
which can not be mimiced by any other foreground, and 
thus might be used in order to devise a method for their de- 
tection. Although most of the current observational efforts 
are dedicated to the study of the (temperature and polariza- 
tion) angular fluctuations, it is clear that measuring these 
spectral features would constitute an unique way to test our 
understanding of the recombination process. For example, it 
will provide a direct determination of the redshift of recom- 
bination, and independent determinations of cosmological 
parameters such as the baryon density. One should mention 
that the broadening due to the scattering off free electrons 
in the epoch of re combination is at most of th e order of 
a few percent (e.g. IDubrovich fc Stolvarovi ri997). implying 
that the features we are discussing in this paper are not 
going to be wiped out. 

Cosmological recombination cannot produce many 
more than one photon per transition in the hydrogen atom 
per one act of recombination. This fact, together with the 
low entropy of the Universe, are the reasons why only in 
two extreme cases might the spectral distortions connected 
with the cosmological hydrogen recombination become ob- 
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servable in the future: (i) the Wien part of the CMB spec- 
trum (Lyman-a, 2s two-photon decay emission & the Balmer 
lines) and (ii) the Rayleigh- Jeans part. 

Concerning the first case, the C1B appears six orders of 
magnitude above the signals we are discussing, so the detec- 
tion in the Wien part of the spectrum will be challenging. On 
the other hand, at very low frequencies the foregrounds like 
the 21 cm line, synchrotron emission of relativistic electrons 
of radio and ordinary galaxies become important and make 
observations more complicated. From this point of view the 
spectral band v > 1.4GHz is more interesting, because the 
21 cm emission will not contaminate in this frequency range, 
and the relative amplitude of the distortions in this band be- 
comes larger (of the order of ~ 10 -7 ). 

The paper is organised as follows. In Sec. 2 we briefly 
review the asumptions of the computations of this paper, as 
well as the basic equations. Sec. 3 presents the results for the 
complete spectrum and for individual lines. In that section 
we discuss the robustness of the results when changing the 
total number of shells, as well as the non-equilibrium effects 
arising from the angular momentum states. Sec. 4 and 5 
present the discussion and the conclusions, respectively. 



2 BASIC EQUATIONS 

In order to calculate the time evolution of the populations 
of the mult i-level hydrogen ato m we adopt the procedure 
described in lSeaeer et alJ i2000l) . Here we would like to em- 
phasize the most important assumptions that this approach 
implies: 

(i) The observed deviations of the CMB s pectrum from 
a blackbody are smaller than AJ/Jo < 10 -4 jMather et alJ 
ll994l:lFixsen et alJll996h . Therefore the radiation field is as- 
sumed to be described by a perfect blackbody with temper- 
ature To = 2.725 K. 

(ii) T he quasi-static solution for th e spectral-line profiles 
is valid dRvbicki fc dellAntoniolll994ft . 

(iii ) The Sobolev escape probability method iSobolevI 
1960) is adopted to deal wit h the evolution of all lines (see 
iDubrovich fc Grachev! (l2004h for a derivation of this method 
in a cosmological context). In practice, this will only af- 
fect those transitions which are directly connecting with the 
ground state, while for the rest the optical depths t v are 
so small that the escape probability is very close to unity. 
Therefore it is possible to work in the optically thin limit 
taking into account only terms proportional to t v for the 
corresponding transition iRubino-Martm et al.l l2005'l. 

In our computat i on, th e major differences with respect to 
the lSeaeer et alJ (|2C)00) procedure are: 

(i) We split the hydrogen levels up according to their 
principal and angular momentum quantum numbers (n, I) 
and include all sub-levels. To illustrate the dependence of 
the recombination spectrum on the treatment of the an- 
gular momentum substates, we also consider some cases, 
in which we follow the populations of all sub-levels with 
n ^ n sp , while for levels with n sp < n ^ n max we aver- 
age over the angular momentum states (i.e we assume that 
the sub-levels are populated according to their statistical 
weights, so that N n i — (21 + 1) iV ns is fulfilled). This implies 



that for given n sp and n max we follow the populations of 
N = n sp [n Bp + l]/2 + 7i max — n ap separate hydrogen states. 

(ii) We neglect collisional rates, as they hav e been found 
to be of minor importance dSeager et al.ll20obTl . 

(iii) We only follow the Hell ionization fraction without 
any detailed follow-up of separate levels. 

(iv) Instead of directly tracking the physical population 
Ni of a level i we define the variables Xi — Ni/Nu, where 
Nn(z) is the total number density of hydrogen nuclei in the 
universe. In this way one can absorb the expansion term. 

A very simila r trea tment has been used in 
iRubino-Martm et alJ (J2005) to calculate the imprints 
of the cosmological hydrogen recombination lines on the 
CMB angular power spectrum including all sub-levels up 
to n max = 10. Here we developed two independent codes 
using different routines of the NAG-Library 2 to solve the 
coupled-system of stiff ordinary differential equations. The 
results of both codes were absolutely consistent and robust. 

In addition, we explored three different numerical 
methods to compute the photoionization and photo- 
recombination rates: the approximations (valid for small 
values o f the energy of the ejected electron) given by 
Burgess] dl95Sl); the poly nomial expressio ns given by 
Karzas fc Latter! il96lTl and lBoanimjm| j|r96j|h an d the nu- 



merical subroutines provided bv IStorevfc^Hummerl il991 
The second and third approaches were found to give identi- 
cal results for all levels up to w max = 30. Even though the 
approximations adopted bv lBurgess! il958T> are not valid for 
high values of n and frequencies far away from the photoion- 
ization threshold of each transition, we found that for the 
calculations of this paper (n max 30), they yield similar re- 
sults (within a few percent accuracy) with respect to those 
obtained in the other two cases. It also turned out that the 
inclusion of induced photorecombination is important for the 
brightness of those transitions with An = 1 at high levels 
(n 3> 1). Finally, we mention that we have also explored 
simple formulae for the computation of the photoioniza- 
tion/photorecombination cross-sections, as those obtained 
using the Kramers approximation, in which the Gaunt fac- 
tors are neglected. We found that this kind of approxima- 
tion is not reproducing properly the trends within a given 
shell, overestimating the rates which are connecting those 
sublevels with very high Z-values to the continuum. 

For all the computations presented in this work the fol- 
lowing values of the cosmological parameters were adopted 
jBennett et alJl2003D : Sl b = 0.0444, Q to t = 1, O m = 0.2678, 
fl A = 0.7322, Y p = 0.24 and h = 0.71. 



2.1 Hydrogen spectral distortions due to radiative 
dipole transitions 

In order to calculate the spectral distortions arising from 
the hydrogen dipole transitions we use a 5-function approx- 
imation for the line-profile, i.e. we neglect any lifetime ef- 
fects (Lorentz-profile) or thermal broadening (Voigt-profile). 
A straightforward derivation yields the spectral distortion, 
/S.Iij(u), arising form the transition of the upper level i to 
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Hydrogen Recombination Spectrum 
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Figure 1. Full hydrogen recombination spectrum for n max = 30. The populations of all the angular momentum substates were taken 
into account. We indicated the contributions corresponding to the Lyman, Balmer, Paschen and Brackett series. The straight line shows 
one possible asymptotic behaviour at low frequencies when including higher hydrogen shells (see text for details). 
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Figure 2. Relative spectral distortion of the CMB spectrum for the hydrogen recombination spectrum as computed in this paper. We 
use the same frequency range as in Fig. Note that the Lyct and two-photon decay features produce very strong relative distortions 
which exceed unity, so they are not presented here for clarity. 
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the lower level j (e.g. also see IWong et alJl2006fr 



ch 



ARij (z e 



(1) 



47T H(z em )[l + z cm ] : 

where z cm is the redshift of emission, which relates the 
observing frequency v to the transition frequency vtj by 
v — + z cm ]. Furthermore H denotes the Hubble- factor 

and the effective radiative rate for the considered transition 
is given by 

A AT. „h.Uij/kT~, 



ARi 



Pa 



l - 



9 j ^ 



-hvij/kT-, 



(2) 



where pg is the S obolev-escape probability as defined in 
ISeager et, alJ l|2f)f)ffi , Aij is the Einstein- ^-coefficient of the 
transition, Ni and (?; are the population and statistical 
weight of the upper and Nj and gj of the lower hydro- 
gen level, respectively. Furthermore we have assumed that 
the ambient photon field is given by a blackbody spectrum 
with temperature T 7 = Tp[ l + z], where Tp = 2.725 K is the 
present CMB temperature jFixsen fc M ather 20q2), and we 
have made use of the Einstein-relations. 




Figure 3. Hydrogen recombination spectrum for different values 
of fimax. The contribution from the 2s two-photon decay was not 
shown on this plot. We show the cases of n max = 10, 15, 20, 25 
and 30, always with n Bp = n max . 



2.2 Hydrogen 2s- Is spectral distortion 

The spectral distortion arising due to the two-photon decay 
of t he hydrogen 2s level can be obtained by the integral (e.g. 
see IZeldovich et alJll968t IWong et al1l2006t) 



A/ 2 slsM 



chv 
4tt 



f 

Jo 



AR 2sls (z)<t>{v[l + z]/i> a ) 
H{z)[l + z] 3 



dz. 



(3) 



Here Ai? 2s is = A 2s i s [N 2!l - N ls e- h " a/kT -'], where the total 
two-photon decay rate for the hydrogen 2s-level is given by 
A 2sls = 8.22458 s" 1 and v a w 2.46 x 10 15 Hz is the Lyman-a 
photon frequency. 

A definition of the profi le function 4>(y) can be found in 
ISpitzer fc Green stein| (| l95lf> but a sufficient app roximation 
was obtained bv lr^jgsbaumer fc Schmutzl jl984li 



4>{y) 



c 



[w(l 



+ a w 



/3+7 Al 



(4) 



where w = y[l - y], C = 24.5561, a = 0.88, j3 = 1.53 and 
7 = 0.8. Here it is important to note that <f>(y) is normalized 
to f Q " a 4>(v/v a ) du — 2 in order to account for the fact that 
two phot ons are emitted per decay o f one hydrogen 2s-level 
(see also lBoschan fc B iltzingcr 1998). 



3 RESULTS 

3.1 The full recombination spectrum 

In Figures Q and |21 we present the full hydrogen recombina- 
tion spectrum for n max = 30. The population dependence 
on I was taken into account for all shells (n sp = 30). In this 
calculation the evolution of 465 separate hydrogen states in 
the redshift range of 500 ^ z ^ 3500 was included. 

At high frequencies one can clearly see the features cor- 
responding to the Lyman, Balmer, Paschen and Brackett 
series. We will discuss some of these lines in more detail be- 
low. At lower frequencies the lines start to overlap strongly 
and eventually merge to a continuum at very low v. The 
figure also shows one possible asymptotic behavior at low 
frequencies. This is meant to illustrate the fact that, when 



including a larger number of shells, the additional cascading 
electrons will enhance the emission at low frequencies. Here 
we decided to plot this asymptotic limit using the slope of 
iKholupenko et all l|2005h . These authors were the first to ob- 
tain the spectral distortions for a very high number of levels 
(n max = 160), although under the assumption of quasistatic 
evolution and without taking into account the populations of 
the angular momentum substates seperately. Their average 
slope in the 1 GHz to 10 GHz range is approximately +0.35, 
while in our case, the average slope in the frequency range 
2 — 10 GHz has a value of approximately +0.8. A detailed 
discussion of this asymptotic region requires the inclusion 
of several physical effects which have not been treated in 
this paper. For example, it is clear that at even lower fre- 
quencies free-free absorption will start to erase any emission 
from the very high hydrogen levels, and in addition colli- 
sions should become important. Therefore, the study of the 
low-frequency asymptotic region will be left for a subsequent 
paper. However, we exp ect that the value o f the slope ob- 
tained from the paper of Kholupe nko et alJ l|2005f) provides 
a lower limit. 

Table Q summarises some parameters (peak amplitude, 
central frequency and relative width) of the first five fea- 
tures in the complete (including the 2s decay line) spectrum 
presented in Fig. Q It is interesting to compare these val- 
ues with those for the corresponding a-transition (An = 1). 
Clearly, the main contribution to the features comes from 
this line, although the emission due to the 2s-decay is sig- 
nificant in the range of the Lyman-series. Furthermore, the 
width of the features is slightly increased with respect to the 
a-line due to the addition of higher transitions (An > 1), 
and in the case of the Lyman and Balmer features, because 
of the 2s decay emission. 



3.1.1 Dependence on the number of shells 

At frequencies v <; 2 GHz the slope of the computed spec- 
trum changes significantly. This is due to the lack of diffu- 
sion of electrons from levels above n max and transitions from 
those levels with An > 1, which have not been accounted 
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Table 1. Positions and amplitudes of the main features shown in Fig. Q We show the central frequency [vq) and wavelength (Ao) as 
observed today, the relative width at half maximum (Av/v) and the peak intensity for the features corresponding to the first five series 
of hydrogen. For comparison, the last five columns show the central redshift of formation, the central frequency, the relative width, the 
peak and relative amplitude of the corresponding a-transition (i.e. An = 1) when considered separately. 
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7.3 x 10~ 28 
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Figure 4. Non-equilibrium effects of the angular momentum 
states on the hydrogen recombination spectrum. We consider the 
case n max = 30 for two different values of n Bp as indicated in 
the legend. The contribution from the 2s two-photon decay was 
omitted here. 



for here. In order to estimate this effect we have performed 
additional calculations with n max = 10, 15, 20 and 25, all 
for n sp = n max . A compilation of the results is shown in 
Figure One can see that the lines of the Lyman, Balmer, 
Paschen and Brackett series basically do not change when 
increasing the number of shells beyond n max ~ 20. On the 
other hand, the low frequencies part (y <C 10 GHz) is still 
varying rather significantly. Our calculations suggest that 
for n max = 30 at the level of a few percent the spectrum 
is converged for v > 20 GHz even when adding more shells. 
Pushing n max to larger values should lead to an increase of 
the distortions at low frequencies and eventually should par- 
tially fill the gap with respect to the asymptotic behaviour 
indicated in Figure 



3.1.2 N on- equilibrium effects of the angular momentum 
states 

Most of the computations of the hydrogen recombina- 
tion spectrum published to date were either based on 
the e f fective three-level atom (e.g. iBoschan fc Biltzingerl 
Il998l : IWong et al.1 l200fj) . or used averaging over angu- 
lar momentum sub-leve ls (e.g. iDubrovich fc Grachevll2004t 
iKholupenko et al"ll2005l) . Here we shall discuss the influence 



of the latter simplification on the solution obtained for the 
recombination spectrum. 

Fig- El compares our result for n ap = 30 with the case of 
n sp = 2 (which represents the standard computation with- 
out detailed follow-up of the Z-substates except for the 2s 
and 2p levels 3 ). The " non-splitting" solution has a higher 
intensity at frequencies smaller than v ~ 2 GHz and also a 
smaller slope in this region. This shows that, for a hydrogen 
atom with given n max , the n ap = 2 solution is effectively 
having a larger number of transitions in the upper shells 
(i.e. higher diffusion), resulting in an excess of emission. We 
will come back to this point below. In addition, there are 
also effects on the high frequency part of the spectrum. For 
n sp — 2 the peak intensities of the main transitions (An = 1) 
are under-estimated while the ratios of the main to the sec- 
ondary transitions (An > 1) are decreased. In other words, 
the incomplete treatment of the high levels also affects the 
diffusion to lower levels, and hence the line intensities are af- 
fected as well. For example, the Balmer lines of the spectrum 
are strongly affected. The ratio of the two peaks appearing 
at ~ 350 GHz and ~ 550 GHz (which correspond to Ha and 
a superposition of lines higher than H7, respectively; see the 
Section on Balmer lines) decreases from roughly 6 to about 
2, and the height of the smaller peak is more than 2 times 
bigger than in the full calculation. 

We note that a co mparison of our res u lts fo r n sp = 2 
with those obtained by IKholupenko et all i2005l) shows a 
very good correspondence. When doing the computations 
using their cosmological parameters, we find good agreement 
also in the amplitudes, showing that our two implementa- 
tions can produce similar results. Given that these authors 
produced their spectrum assuming a quasistatic-evolution 
for all shells above n = 2, our results further supports the 
idea that their approach is appropriate (see also the com- 
ments in IDubrovich fc Shakhvorostoval 2004). However, in- 
cluding the evolution of the populations of all Z-substates in 
this approach seems more difficult to solve, because the de- 
viations from statistical equilibrium within a shell are very 
small, and the corresponding columns in the matrix of the 
level populations will be degenerate. To test this point, we 
have also implemented the quasi-static approach for shells 
n > 2, while solving the evolution (differential equations) for 
n ^ 2 and the electron fraction. For our implementation of 

3 Hereafter, we will refer to the case n sp = 2 as the "non- 
splitting" solution, and to the case n Bp = n m a X as the full follow- 
up case. 
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Figure 5. Non-equilibrium effects on the populations of the angular momentum substates for a given n. The results were obtained 
for n max = 30 following the populations of all the angular momentum sub-states. We present the ratio N n i/Nf^ , where the statistical 
equilibrium (SE) population is computed from the actual total population of the shell by = [(21 + l)/ra 2 ]iVtot- Six cases are shown, 

corresponding to values of the principal quantum number n = 3, 5, 10, 15, 20 and 25. 



the problem, we found that, when considering all the angu- 
lar momentum substates separately, with standard methods 
solving the complete set of ordinary differential equations 
seems more efficient (in terms of computational time) than 
solving the corresponding algebraic system. 

Finally, we have also performed several computations 
considering values of n sp in the range between 2 and n max , 
mainly with the aim of understanding whether the full 
computations could be simplified. In these cases, we found 



that the spectrum obtained is closer to the full solution 
(n sp = n max ) in the high frequency part. However, there 
always appears an abrupt change of slope in the spectrum 
precisely at the frequency corresponding to the main tran- 
sition of the n ap -series. Below this frequency, the average 
slope of the spectrum is natter than in the full solution. 

Obviously, all the effects discussed in Fig. 0] appear be- 
cause the inclusion of i-substates modifies the populations of 
the levels. In order to understand these modifications, and 
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how the differences in the populations of the levels are con- 
nected with changes in the spectral distortions, we shall now 
discuss these effects in detail. 

Effects within a given shell. We first discuss the situa- 
tion for a given shell (i.e. we fix n), where we compare the 
populations of the different sub-levels with the correspond- 
ing expected value in the case of statistical equilibrium (SE) 
within the shell, i.e. when JV^f = [(2/ + l)/n 2 ]iVtot, where 
iVtot is the (actual) total population of the shell. 

In Fig. |3 we present the ratios N n i/N„f of the true 
population N n i of the level for the full computation with 
Wmax = 30, to the expected SE population. We show six dif- 
ferent cases, corresponding to values of the principal quan- 
tum number n = 3, 5, 10, 15, 20 and 25. For each one of 
these cases, we present the ratios at three different redshifts 
(z — 1200, 1300 and 1400) which correspond to the region 
where the lines are forming. At higher redshifts, the devia- 
tions from statistical equilibrium show a similar trend, but 
they are much smaller in amplitude and tend to dissappear. 
This behaviour is expected since one is then approaching the 
epoch of full thermodynamic equilibrium iZeldovich et alJ 
1968). In our calculations the non-equilibrium effects are 
compatible with zero within the numerical precision of the 
two codes at z J> 2000. 

The general trend for all shells above n = 5 can be 
summarized as follows: within a given shell the sub-levels 
with I < 5 — 6 are underpopulated as compared to N^f, 
with the difference being the strongest for the d-level (I = 2) 
and smallest for the level I ~ 5 — 6. This shows that in 
the calculations with full follow-up these angular momentum 
states depopulate faster. For I > 6 the population of the sub- 
states is larger than in the SE situations 4 . This implies that 
these angular momentum states depopulate slower, which 
is likely due to blocking of the transition to lower levels 
(Al = ±1 restriction). On the other hand, for those shells 
n ^ 5, the p and d-levels are also underpopulated, but now 
the strongest deviation appears for p-level (I = 1). When 
comparing our results for n max = 30 with the case of n max = 
25 and n max = 20, we can conclude that the results in Fig. [S] 
for n ^ 10 have already converged, but higher shells could 
increase their deviations from SE, especially at high i-states. 

Although the explanation of these non-equilibrium ef- 
fect requires one to study the net rates for all transitions 
connecting to that particular level, we can understand some 
of the general trends with the following arguments. It is 
clear that assuming SE is equivalent to "instantaneous" re- 
distribution within a given shell, i.e. if an electron is cap- 
tured by a given Z-state, this will immediately increase the 
population of all the other Z-states. One would expect that 
instantaneous redistribution will operate when collisions 
(n, I) — > (n, I') are effective. This should eventually happen 
in sufficiently high shells (see Sect.^J. On the other hand, 
when instantaneous redistribution is not working, then the 
depopulation of each level will depend on the exact route 
the electron can take to reach the ground state. 

As the Universe expands, the blackbody radiation field 



4 Note that this is not true for the very high i-substates and the 
outer shells, as n = 20 or n = 25, as a consequence of the lack of 
cascading electrons from higher levels. 
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Figure 6. Non-equilibrium effects in the populations of the 2s 
and 2p levels. We present the ratio of the 2s and 2p populations 
referred to the equilibrium ratio (i.e. the statistical weights), as 
a function of redshift. Although the deviations are very small, 
they indicate that there is a deficit of electrons in the 2s state as 
compared to the 2p state as a consequence of the recombination 
dynamics. 



is redshifting to lower energies, and thus the number den- 
sity of high energy photons decreases with time. When the 
radiation field does not contains enough photons to main- 
tain equilibrium in a given shell n (i.e. when the number 
of photons with energy above the photoionization potential 
of the shell \n falls below 1 per atom), then the photoion- 
ization rates become smaller than the photorecombination 
rates and the electrons can not escape to the continuum. 
Note that this is first happening for lower shells with the 
largest photoionization potential. 

In the above situation, the electrons try to cascade down 
to lower levels. Our computation only includes radiative 
rates, so electrons in a given shell can only diffuse down 
through transitions with Al = 1. In general, it seems that 
the fastest transitions are found to be those which can di- 
rectly connect with the 2s and 2p-levels, so it is reasonable 
to expect that all sub-levels with I < 2 will be able to depop- 
ulate faster in the complete (n sp = n max ) computation, and 
their SE populations will seem larger. On the other hand, 
due to the bottle-neck, which is produced in the lower states 
during cosmological recombination, states with high values 
of I depopulate slower in the full computation, since elec- 
trons have to cascade down via many intermediate shells. 

Finally, we also present in Fig.HJthe ratio of the popu- 
lations of the levels 2s and 2p as a function of redshift. This 
figure quantifies the fact that at high redshifts, both the 2s 
and 2p populations are in statistical equilibrium, so their 
population ratio equals the ratio of the statistical weights. 
However, as the recombination proceeds, deviations from the 
equilibrium appear, in the sense that there is an overpopu- 
lation of the 2p level with respect to the 2s level because the 
main channel of recombination at low redshifts (z < 1400) 
is the 2s two photon decay. We will return to this point in 
Sec 3.2 and 3.3. 

Effects on the total population of different shells. 

The effects connected with the departure from SE within a 
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Figure 7. Non-equilibrium effects in the populations of the levels 
for different values of n. We present the relative change in the total 
population of the levels (TVtot (n) = ; N n i ) due to the detailed 
follow-up of the angular momentum sub-states, for a hydrogen 
atom with n max = 30. Note that the fractional change, defined as 
r (n S p=nmax) /w-( n =p= 2 )l /i\ r ( nB P= 2 ) 



IK 



'l/K 



is always positive, so that 



the total population of a given shell in the calculations for n sp 
n max is always larger than for n Bp = 2. 



given shell are not only visible for different sub-states, but 
they also produce a net effect on the overall population of 
the level. We illustrate this point in Fig.Q where we present 
the relative change of the overall population of a given shell 
for n max = n ap , with respect to the the standard (n sp = 2) 
computation. 

For a given shell n, the total population (sum over all I 



states) in the computations with n s 



is always larger 



than for n sp = 2. This shows that in general the overall 
depopulation of a given shell is slower for full follow-up. 
Moreover, for n > 3 the relative difference becomes larger 
when increasing the value of n (for n = 3 the net effect is 
slightly smaller than for level n = 2). For high values of n 
(close to n max ), this relative change seems to decrease again 
as a consequence of the lack of diffusion from higher shells. 

At low redshifts z < 900 the relative difference becomes 
larger than 1% for all the levels. On the other hand, in the 
redshift range where the spectral distortions are produced 
(e.g. for the main transitions, the redshift of formation is 
2 ~ f300 for high n, and z ~ f400 for n = 2), the frac- 
tional change is extremely small. This shows that the total 
population of the shell during the epoch of recombination 
does not depend significantly on the treatment of the angu- 
lar momentum sub-states, although these small deviations 
are sufficient to produce the effects shown in Fig. 2] 

Effects on the transitions between different shells fn 

this paragraph we wish to understand the effects on the tran- 
sition rates between different shells. Although the detailed 
analysis again would require one to show the differences in 
the net rate for each transition, we will show that using the 
results and arguments of the last two paragraphs we can 
understand the general trend and sign of the correction to 
the resulting spectral distortion as shown in Fig. 0] 

When including /-substates, the populations of all levels 
will be given by AT, = N? s + ANi, where we introduced JVf s 



referring to the populations in the case of 'non-splitting' in 
angular momentum above n = 2, i.e. n sp = 2. Given that the 
optical depth for al l transitions to levels above the ground 
state is very small jRubino-Martm et al.ll2005l) . the escape 
probabilities pij are very close to unity and therefore will 
practically not change when following the populations of the 
Z-substates. Thus, the net rate can be written as ARij = 
ARfj [1 + A], where ARfj is the rate as obtained in the 
calculation for n sp = 2. The function A is defined by 



ARf? x A 



AiVj 
AT 



(ANj/N? 
(ANj/N? 



Si N? 



-hvij /kT-y 



(5) 

For simplicity, let us now only consider the a-transitions. 
These give the main contribution to the CMB spectral dis- 
tortion and always lead to positive features (see Fig^J, im- 
plying that 



ARij <x 



i _ ii 



-hVij/kTy 



> 0. 



(6) 



Therefore the sign of A will be defined by the rhs of Eq. @. 

Now, comparing Eq. © with (|SJ one can conclude that 
for (ANi /Ni 3 ) /(ANj /N™) > f the term within the brack- 
ets is positive, and becomes negative if this ratio is smaller 
than f. fn addition one has to check the sign of the factor 
ANj/Nf s to obtain the overall sign of A. 

A detailed analysis would require one to study the sign 
and amplitude of these corrections for every one of the dif- 
ferent channels contributing to a given transition. These 
numbers can in principle be derived by combining the in- 
formation of Fig. |21 and [T] However, it is clear that a better 
approach would be to directly study the net rates for ev- 
ery transition of interest. We will do that in the following 
subsections for some examples, although here we want to 
illustrate that one can understand the overall sign of the 
correction to the spectral distortions with this argument. 

If we apply Eq. |^| to the total populations of the two 
levels that we are connecting (this means that we neglect the 
deviations from SE within the shell as have been discussed 
above) , then we can use Fig. [7] to understand the net sign 
of the correction. All the fractional changes in Fig. Q arc 
positive, so we only need to consider the ratio between the 
levels involved in the transition. On one hand, for 2 < n < 20 
the relative population is growing with n, so Eq. |S] shows 
that the correction term to the net rate should be positive, 
and thus the spectral distortion in the full computation lies 
above the result for non-splitting. However, for high n-values 
(close to n max ), Fig. □ shows that (A^/N™) /(AN 3 /N™) 
can be expected to be less than unity for n > 20, implying 
that the term in Eq. |S] inside the brackets can eventually 
become negative, and thus the spectral distortion in the full 
computation will lie below the non-splitting solution at low 
frequencies. 

We note that the above argument is only used here to 
show the qualitative behavior of the correction when tak- 
ing into account angular momentum sub-states. In our real 
computations, we considered all possible radiative channels 
(and hence all possible sub-levels), which contribute to a 
given transition, so a detailed analysis based on the net rates 
ARij for all channels would in principle be possible. 
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Figure 8. The Lyman series and the the 2s two-photon decay 
spectrum. The main contribution to the Lyman series (thick solid 
line) comes from the Lyman-o-transition (dashed line). The peak 
intensities for the first few transitions are AI V = 4.8x 10 — 27 , 2.9x 
1CT 29 , 4.9 X 1CT 30 and 2.2 X 1CT 30 Jm" 2 s" 1 Hz" 1 sr" 1 for the 
Ly-ct, Ly-/3, Ly-7 and Ly-5, respectively. For comparison, the 
2.725 K blackbody is also shown (thin solid line). 



Figure 9. Redshift of formation for different lines. We present 
the normalised intensity for several transitions with An = 1 as a 
function of redshift 2. For comparison, we also show the Thom- 
son visibility function normalised to unity at the peak, showing 
that the redshift of formation of the lines is higher than that of 
formation of the CMB fluctuations. 



3.2 The Lyman-a line 

As found by other authors, the Ly-a line leads to the 
strongest distortion of the CMB spectrum and appears in 
the Wien tail of the blackbody spectrum. In our compu- 
tations it peaks at wavelength around 170 /jm (see Fig. |H| 
and Table and is very similar in shape a nd amplitude 
to the results obtained bv lWong et, al.l l)200(tl . However, as 
discussed below, we do not find any pre-recombination peak. 

In Fig. [5] we present several transitions with An = 1 
in redshift space. The figure shows that the bulk of the 
Ly-a emission origin ates around z ~ 1400, as obtained 
by other authors (e.g. iLiubarskii fc Sunvaevll983l) . Increas- 
ing n one can see that the corresponding transition, due 
to the decrease in the transition energy, stays longer in 
equilibrium with the photon field. However, when the tran- 
sition frequency close to the decoupling redshift (i.e. the 
redshift at which the corresponding line starts appearing) 
drops below the maximum of the CMB spectrum (here in 
terms of photon number, N v oc v 2 /[exp(/i^/fcT 7 ) — 1]), i.e. 
v m ax ~ 90.3 x (1+z) GHz, then increasing n further will shift 
the redshift of formation for the corresponding line back to 
larger z (see Fig. ^transition 10 — > 9). For a more quantita- 
tive estimate one also has to include the dependence of Aij 
on n and eventually will end up looking at the net transition 
rates in more detail. 

Here, another important point is that observations of 
these lines permit us to explore an epoch of the Uni- 
verse, where the optical depth due to Thomson scat- 
tering is very high. Note that the effects discussed in 
iRubino-Martm et al.l i2005f) also provide spectral informa- 
tion but, due to the shape of t he Thomson visibility function 
JSunvaev fc Zeldovich|[l970al) . this information comes from 
lower redshifts (z ~ 1100). 

Fig. 1101 shows the redshift dependence of the net (ra- 
diative) rates of the Ly-a and 2s two-photon decay transi- 
tions. As it is well-known, the 2s— Is rate dominates during 
the epoch of recombination, and provides the main channel 
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Figure 10. Net rates for transitions from levels n = 2 to the 
ground state. The rate AiJ(2p — Is) mimics the shape of the Ly- 
a line (compare with Fig.[5J. Note that the 2s two- photon decay 
rate dominates for redshifts smaller than ~ 1400 (the two rates 
are equal at z ~ 1380), i.e. where the bulk of recombination takes 
place, and where all the contributions from higher transitions 
originate. 



throu gh which recombination can proceed iZeldovich et alJ 
1968). On the other hand, escape through the Ly-a line hap- 
pens at earlier times (the two rates are equal at z ~ 1380). 
Note that due to the narrowness of the Ly-a line-profile, the 
shape of Ai?(2p — Is) mimics that of the Ly-a contribution 
to the CMB spectral distortion (compare with Fig. |UJ . 

3.2.1 About the pre-recombination peak 

In earlier calculations of the recombination spectrum the 
existence of a pre -recombination peak for the Lyman-q tran- 
sition was found jRvbicki fc deU'Antoniclll99'l IWong et alJ 
2006). However, in our calculations we do not find any sig- 
nificant emission from any hydrogen line above redshifts 
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z ~ 2000, provided that the primordial radiation field is 
described by a pure blackbody law. 

With the assumption that there are no intrinsic dis- 
tortions of the CMB spectrum it is expected that at some 
high redshift the populations of all neutral hydrogen states 
can be described using the Saha-relations, with tiny devia- 
tions only coming from the disturbance of equilibrium due 
to the expansion of the Universe. Clearly, in this case it is 
impossible to have any uncompensated loops between dif- 
ferent levels and hence no signific ant net emission should 
occur llLiubarskii fc Sunvaevl ll983'l, Our calculations show 
that, at redshifts z > 2000, the net rate of the 2p— Is process 
becomes zero (within the numerical accuracy of the codes, 
which is much higher than the effects we are discussing here) . 

We found that here it is very important to obtain the 
solution for the populations of levels with sufficiently high 
accuracy. In local thermodynamic equilibrium (LTE) the 
Saha-relation yields J 



_ 9j_ e h Vij /kT^ 



(7) 



for the ratio of the population of the levels i and j. This 
ratio evidently is very sensitive to the value of the transition 
frequency Vij . Inserting this expression into Eq. it is clear 
that analytically one obtains no net emission, provided that 
v ij — \Xi ~ Xj\/hi with Xi being the ionization potential of 
the corresponding level. 

However, in numerical calculations at sufficiently early 
epochs one will obtain a solution for the population 
of each level, which is extremely close to the numeri- 
cal LTE solution. But in fact, when one now computes 
fr(N j /N i )l™e- hvi i /kT ~< numerically even deviations of 

the order of A (Nj /Ni)^™ / {N 3 /Ni) LTE ~ 10~ 5 from the 
real LTE solution will lead to sizeable false emission, when 
deviations from LTE should be small. 

This in principle can occur when using a slightly differ- 
ent value, Vij — + Svij (for example from older precal- 
culated tables or using different combinations of the natural 
constants for unit conversions) for the transition frequency 
in the setup of the coupled system of ordinary differential 
equations, while obtaining the output of the recombination 
spectrum with . In other words, a false pre-recombination 
peak appears when one is not comparing the obtained nu- 
merical solution with the corresponding numerical LTE so- 
lution, which in this setup should fulfill 



N, 



9j Mvij+Sv^/kT^ 



Ni I Qi 
1 / num 3 L 

Inserting this expression into Eq. one obtains 
&mf m {v)<xl-e hS ^ /kT -< . 



(8) 



(9) 



Because of the large factor k — hvij/kTo in the expo- 
nent even for Svij/vij ~ 10 -6 this leads to notable false 
emission at early epochs, while leaving the main recombi- 
nation peak unchanged. As an example, for the Lyman-a 
transition (k ~ 43455 and Av a /v a ~ 10 -6 ) one obtains a 
pre-recombination emission, which peaks at z ~ 3100 with 
AI ~ 6 x 10" 26 Jm -2 s" 1 Hz" 1 sr" 1 , while the recombina- 
tion peak is roughly 10 times smaller! 



Note that in local thermodynamic equilibrium T c 



One should mention that although at z > 2000 there is 
a small difference in the electron and photon temperature 
(AT/T-y ~ 10 -6 ) this cannot lead to any pre-recombination 
emission: due to the large entropy of the Universe (Nh fN 1 ~ 
10~ 10 ) one expects that even for a larger AT/T 7 the statis- 
tics of the hydrogen levels is controlled by the photon field. 
This implies that in the Boltzmann factors T 7 appears in- 
stead of T c . It is expected that even the inclusion of col- 
lisional processes, which we neglected here, will not change 
this fact. Similarly, the expected small change in the electron 
temperature due to Compton heating a distorted photon 
field will not be of importance. We therefore conclude that 
there should be no hydrogen pre-recombination emission and 
that similarly one should not expect any pre-recombination 
emission from helium. 



3.3 The two-photon decay spectrum 

In Figure |H| we also show the contribution due to the 2s 
two-photon decay. The peak of the emission is roughly at 
v ~ 1460 GHz (w 205 /im) and partially fills the gap between 
the Lyman and Balmer series (compare with Fig. 0. The 
amplitude of the 2s two-photon decay line is approximately 
3.3 x 10~ 27 J m~ 2 s _1 Hz -1 sr -1 , which is only a factor ~ 1.5 
times smaller than the peak value of t he Ly-a line. This ratio 
is com patible with the result found bv lBoschan fc Biltzinger 
lll99Sf) . while it differs from the result of IWong et al.l d2006ft 
by a factor of 2, as a consequence of their different normal- 
isation of the profile function. 

Fig. 1101 also presents the 2s— Is net rate as a function 
of redshift. Note that one can not directly infer from this 
figure the position of the peak intensity of the contribution 
to the CMB spectral distortion due to the 2s decay, because 
in this case the line profile (which appears in the kernel of 
integration in Eq.0 is very broad. Therefore, the resulting 
spectral distortion is much broader than any other radiative 
line (Av/v » 0.8, which is practically four times larger than 
the value for Lya), and it peaks at higher frequencies than in 
the simple ex trapolation of the redshifted central frequency 
v 2s ~ v hya /2 JWong et alJl2006li . 



3.4 The Balmer transitions and beyond 

Fig. llll shows the different contributions to the Balmer series 
using a linear scale in intensity. The main contribution comes 
from the Ha line. One can see that the contribution from H/3 
becomes negative in a certain frequency range. The positive 
wing of H/3 slightly enhances the high-frequency wing of the 
main feature at 350 GHz, while the rest of the lines (H7 and 
above) add up to produce a secondary peak around 600 GHz, 
completely compensating the negative part of H/3. Note that 
this re sult differs from that of iDubrovich fc Shakhvorostoval 
( 2004), where their H/3 and H7 are both negative, and all the 
higher transitions do not compensate them. Also in general 
the amplitude of our results is sma ller by a factor of ~ 10. 
We al so mention that the results of iRvbicki fc dell' Antonio! 
( 1993) for n max = 10 contain negative features which are 
not found in our computations (see Fig. [HJ. 

Fig. H2l shows the net (radiative) rates for all transitions 
connecting levels n — 3 and n — 2, and which contribute to 
the Ha line. The fastest rate is associated with the 3p— 2s 
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Figure 13. Different contributions to a high-n series. Left panel: we show the contribution of all transitions to level n = 15 for a 
hydrogen atom with n max = 30. There is a strong peak produced by the corresponding An = 1 transition, while the rest of transitions 
(An 2) add-up to produce a (practically) featureless high-frequency wing of the series. Right panel: we show the net contribution of 
all transitions to levels n = 14, 15 and 16, in normalised intensity units. At frequencies below a few GHz, the contribution of the An = 1 
transitions start to overlap. 
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Figure 11. Different contributions to the Balmer series. For con- 
venience we normalised each of the contributions to unity at their 
maximum. The normalisation factors are A/„ = 6.5xl0 -27 , 7.5 X 
1CT 28 , 4.1 X 1CT 28 and 2.5 X 10~ 28 J m -2 s _1 Hz -1 sr _1 for the 
Ha, H/3, H7 and H<5, respectively. 



channel, as one would expect because the Ha line forms 
around z ~ 1300 (see Fig. |5J|, where escape through the 
2s— Is channel is faster than through the 2p— Is one. 

For higher series, one can perform a similar analysis. 
For instance, in the left panel of Fig. 1131 we show the first 
contributions (An ^ 5) to the series corresponding to level 
n = 15. It is clear that for large values of n, the relative 
distance between transitions with An = 1 and An = 2 in- 
creases (Av/v = [fn,n+l — Vn,n+2]/Vn,n+-L > 0.8 for n > 13). 
Given that the typical width of the lines is of the order of 
0.3 (see Table0, the a-transition is giving a separated peak, 
while the other transitions of the series overlap and result 
in a nearly constant emission plateau in the high-frequency 
wing. Although it is not shown here, we have checked that 
the slope of this plateau is steeper than the one obtained in 
the case of n sp = 2. The co-added contribution of all these 



Figure 12. Net rates for transitions from levels n = 3 to levels 
n = 2. This figure contains all transitions contributing to the Ha 
line, namely Ai?(3s - 2p), A_R(3p - 2s) and A_R(3d - 2p), except 
the 3p— Is. However, the contribution of the rate Ai?(3p — Is) 
to the line is very small (it peaks at z ~ 1470 with amplitude 
4.6 X 10 — 7 m -3 s" 1 ), so we can omit it. During the epoch of 
formation of the line, the fastest transition is 3p— 2s, and therefore 
the population of the 3p level will be smaller than the expected 
equilibrium value (see Fig. l5l. 

high-frequency wings of all series is contributing to the net 
slope at low frequencies in Fig.Q 

On the other hand, the peaks coming from the main 
(An = 1) transition of a given series cannot be sepa- 
rated from those of other series for high values of n, be- 
cause they start to overlap. In the right panel of Fig. 1131 
we show this overlapping for the three series correspond- 
ing to n = 14, 15 and 16. In this case, Avjv defined as 
[vn,n+i — Vn+i,n+2\/vn,n+i is a decreasing function with n, 
so Av/v ^ 0.2 for n ^ 13, and thus the peaks start to over- 
lap and merge into a continuum. This effect can be seen in 
the low- frequency part of Fig.0 Note that the co-added con- 
tributions of all the high-frequency wings is giving the net 
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slope. If the populations of Z-substates are not followed, the 
transitions with An > 1 will be larger, and thus we will have 
a different behavior in the asymptotic slope in the sub- GHz 
range. 



3.5 Non-equilibrium effects of the angular 

momentum states on the electron fraction 

One of the important variables in our calculations is the 
electron fraction, x c — N £ /Nh, as a function of redshift z. 
Here we want to quantify the effects on the electron frac- 
tion. To do so, we compare the results obtained for the 25- 
shell atom in the full calculation (n sp = 25) with those for 
standard computation (n sp = 2). Fig. 1141 shows the relative 
change of the free electron fraction due to the inclusion of 
angular momentum substates. In the standard computation 
the residual electron fraction at low redshifts (z — 800) is 
smaller by ~ 1.5% than in the case of full /-follow-up, imply- 
ing that in this case the process of recombination is slightly 
slower. Note that for n sp — 10 one obtains approximately 
half of the effect. 

The net eff ect during the epo ch of recombination 
(zdcc = 1089 ±1, iBennett et afll2003ft is smaller than 0.1%. 
Therefore one would only expect a small impact on the 
angular power spectrum of the CMB. To quantify this ef- 
fect w e have modified the Cmbfast iSeliak fc Zaldarriaeal 
code, to incorporate the recombination history ob- 
tained in our computations. We note that this recombina- 
tion history is comp uted inside Cmbfast using the Recfast 
iSeaeer et aljfl999T) code, which is based on calculations for 
a 300-level hydrogen atom in which SE among the I sub- 
levels was assumed iSeager et aljEoOoh . Our computation 
makes the assumption that the relative change of the elec- 
tron fraction, as obtained for the 25-shell atom, will not vary 
significantly when including more shells. Although the inclu- 
sion of more shells in our calculations should still lead to a 
decrease of the residual electron fraction at low redshifts 
(one would expect to converge for n max ~ 50), we found 
that the above assumption is justified reasonably well. We 
will further discuss this point in a future work. 

Fig. 1151 shows that the relative difference of the Ce 
is smaller than ~ 0.5% (~ 1.1%) for the TT (TE) 
angular power spectrum. Following Sea ger et alJ (Il999l) 
we can model the results by multiplying the recombina- 
tion/ionization rates by a 'fudge factor'. Our results for x e 
are reproduced with an extra factor 0.967, i.e. the fudge 
factor inside Recfast should be 1.10 instead of 1.14. 

Finally, we would like to mention that there are other 
effects, which may produce small c hanges in the Ct's. For 
example, Chlu ba fc Sunvaevl a2006l) considered the effect of 
stimulated emission due to the presence of sof t CMB pho- 
tons o n the total 2s-decay rate. In addition, iLeung et alJ 
(2004) discussed the effect of a modification of the adiabatic 
index in t he matter equation of state a s recombination pro- 
ceeds and lDubrovich fc Grachev! ^200511 included corrections 
due to the two-photon decay of higher hydrogen levels and 
helium. All the aforementioned processes lead to changes of 
the order of 1% and even cancel eachother partially, leaving 
a lot of space for future work and discussion. 



7?, 5p = 7J, mQ , 

n,= 10 



Figure 14. Relative change in the free electron fraction (x e = 
Nb/Nh)- We consider the 25-shell atom and show the rela- 
tive difference in x e with respect to the standard computation 
(Ax c /x c = [x c — x c (n Bp = 2)]/x e (n Bp = 2)) for two cases: full 
computation (n sp = n max ) and n Bp = 10. 




1000 



Figure 15. Relative change in the temperature (TT, solid line) 
and polarization (TE, dashed line) angular power spectra due 
to the inclusion of detailed follow-up of the angular momentum 
sub-states, for the WMAP best-fit cosmological model. 



4 DISCUSSION 

The spectrum of the background radiation of the Uni- 
verse has a minimum i n the vicinity of 30 — 40 cm 
llLongair fc Sunvaevl Il972l) . It is very impressive that the 
brightness temperature in the recombinational line compo- 
nent has a maximum relative to the background radiation 
spectrum precisely in this spectral band (see Fig. |5J, and 
reaches values of the order of 10~ 7 . This is a consequence of 
the fact that the asymptotic slope of the spectral distortions 
to the blackbody (which in our case is of the order of « +0.8) 
is decaying slower with frequency than the Ray leigh- Jeans 
part of the CMB (slope +2). Note that the turn-over at 
around 0.5 GHz is due to the lack of diffusion from high 
levels in our calculation with n max = 30. 

It is extremely difficult to separate an additional con- 
tinuum component with quasi-flat spectrum from the fore- 
ground sources. However, the lines with An = 1 will still in- 
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fluence the shape of the integral curve, and the contribution 
from recombination is modulated in such a way that no fore- 
ground source is able to mimic it. Obviously, one of the most 
interesting spectral regions corresponds to those frequencies 
just above 1420 MHz (21 cm line), where the contribution of 
synchrotron emissio n of relativistic electrons of radio and or - 
dinary galaxies fe.g. lLongair fc Sunvaevll972HLongairll993l) 
and free-free emission from our Galaxy are not expected to 
have a quasi-periodic dependence on the wavelength (for a 
mo re detailed d i scussi o n of foregrounds for 21 cm studies, 
see lOh fc Mackl j2003MMorales et alJ d2005l) ). 

At low frequencies, the amplitude of this frequency 
modulation is not strong, but it might become measur - 
able by future experiments (see e.g. lFixsen fc MatheJ 2002'l. 
Fig. 1161 shows that the contrast of these features in the fre- 
quency range between 1 and 10 GHz is between 5% and 
30%. Even though the amplitude of these features may still 
increase when considering larger values of n max , the con- 
trast should not change significantly. Note that the Balmer, 
Paschen and Brackett series have a higher contrast and 
they are found in the frequency band where the contami- 
nation is smaller, but they will be much more difficult to 
observe because they lie two orders of magnitude below in 
relative intensity (see Fig. |5J. We would like to mention 
here that there have been observational effort s to measure 
the CMB spectrum at ce ntimeter wavelengths l|Kogut et alJ 
120041: iFixsen et al.ll2004l) . searching for signatures of spec- 
tral distortions from early energy injection in the spectrum 
JSunvaev fc Zeldovichlll 97014 iTllarionov fc Svunva.evllT975h. 

Although the Lyman-a distortion is the strongest fea- 
ture in the spectrum (relative distortion in the peak is of 
order unity), the Cosmic Infrared Background (CIB) in this 
region is dominant, making its observation extremely diffi- 
cult. However, the sensitivity of the detectors is improving a 
lot, so these observations could become feasible in the future. 
For e xample, the recent measurements of Spitzer ijDole et alJ 
2006) have shown that at 160 /im, practically three quarters 
of the Mid-Infrared Background (MIB) can be already re- 
solved into luminous infrared galaxies at z ~ 1. If the MIB 
could be completely resolved, one might be able to remove 
the source contamination in a given piece of sky in order to 
measure the spectral distortions discussed here. 

For illustration, Fig. ll7l shows the amplitude of the spec- 
tral distortions compared to the d etermination of the CIB 
obtained wit h COBE FIRAS data (IFixsen et alJll998f) . and 
the result of lLagache et all J2000) in the Lockman Hole us- 
ing the WHAM Ha survey and the Leiden/Dwingeloo HI 
data. We note that, in addition to the CIB signal, there will 
be also a st rong contaminatio n from our Galaxy (see the 
discussion in lWong et al.ll2006T) . 

We would like to stress that all the results of this work 
were obtained within the standard framework which has 
been adopted by many other authors working on this topic 
(see the references in the introduction). The main difference 
with respect to these earlier works is the detailed treatment 
of the i-substates. However, it is clear that one has to eval- 
uate the impact of other physics which has been neglected 
so far. As an illustration, we will discuss three additional 
physical processes which have not been included here and 
which could modify the results. 

The first effect is connected with collisions. Although 
the collisional rates are smaller than the radiative rates dur- 



ing the epoch of recombination ISeager et al.ll200oT) . they 
could produce a redistribution of electrons if we consider 
states which are not directly connected via dipolar electric 
transitions (e.g. collisions can produce transitions of the type 
(n, I) — > (n, I') also for Al 7^ 1). In principle, this effect 
should be important for shells with large principal quan- 
tum number n, and would tend to restore SE, thus restoring 
the (21 + l)-ratios of the populations of the levels within a 
given shell. In other words, one should eventually recover 
the non-splitting solution for sufficientl y high shells, and 
the as ymptotic slope of +0.35 found bv iKholupenko et alJ 
(2005) might be reached. One can do a simple estimate of 
th e effect by using the coll isional rates described by Eq. 43 
of IPengellv fc Seatonl lll964f) . In our problem, the collisional 
transition probability q n iN n iN e becomes of the order of the 
radiative transition probability A n i for n > 25 — 30 at the 
redshift of formation of the lines (z ~ 1300). Thus, if one 
wants to include a higher number of shells in the computa- 
tions, the effect of collisions on the mixing of the populations 
in different Z-sublevels has to be taken into account. 

The second effect is the free-free absorption (see e.g. 
Illari onov fc Svunvaevl Il975t iDubrovich fc Shakhvorostoval 
l2004f) , which at low frequencies will tend to restore the black- 
body with T 7 = T c and thus will erase the features in the 
spectrum. However, this effect will be only important for 
high shells, so our computation up to n max = 30 will not 
change significantly. 

Finally, also the inclusion of additional radiative chan- 
nels (stimulated two-photon decay from higher levels, forbid- 
den transitions, i.e. both electric and magnetic quadrupole 
transitions, or magnetic dipole transitions) may have some 
impact on the computed spectral distortions. However, we 
would expect that the inclusion of these new channels will 
not subtantially modify the spectral distortions discussed 
here. The reason for this is that all the sub-levels with n > 2 
can connect to lower states via permitted radiative transi- 
tions having escape probability very close to unity. There- 
fore, the inclusion of additional "slow channels" should only 
modify our results at the level of a few percent. We have al- 
ready c hecked that the inclusion of stimulated 2s two-photon 
decay dChluba fc Sunvaevll2006T) and some additional two- 
photon decay channel s for higher hydrogenic level s (follow- 
ing the prescription in lDubrovich fc Grachevl2005l) . leads to 
differences with respect to our computations of the order of 
1%, in both the recombination spectrum and the electron 
fraction. 

Certainly one should try to develope an approach in- 
cluding all the aforementioned processes stimultaneously, 
but this shall be left for a future work. 



5 CONCLUSIONS 

In this paper, we have computed the spectral distortions of 
the CMB frequency spectrum originating from the epoch of 
cosmological hydrogen recombination, by following the evo- 
lution of the populations of the levels of a hydrogen atom 
with n max = 30 in the redshift range from z — 3500 down 
to z = 500. Our computation follows in detail the evolu- 
tion of all angular momentum sub-states. This permits us 
to describe the non-equilibrium effects which appear dur- 
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Figure 16. A detail of the hydrogen spectrum in the frequency 
domain between 1 GHz and 10 GHz. We present in linear scale 
the amplitude of the features in the region of the spectrum were 
we expect that detection should be more feasible. Dotted lines 
are connecting the maxima and minima of the features, so one 
can read directly their contrast. 
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Figure 17. Spectral distortions and foregrounds. We compare 
the result of Fig. with the pure blackbody spectrum of 2.725 K 
(dotted line); with the emission of the Cosmic Infrared Back- 
ground derive d bv [ Fixsen et all ll99STl (dot-dashed line) and 
lLagache et alj i2000l) (solid line between 200^tm and 1.2mm); 
and with the estimates of th e contribution of cxtragalactic radio 
sources from iLongairl il993t) (note that this refers to the extra- 
galactic sources with an average spectral index of —0.8 in flux, 
and which roughly scales as 6 K x(v/178 MHz) -2 8 ; this is why 
in these coordinates, their contribution seems practically flat). 

ing the process of recombination in the populations of the 
sub- levels. 

The populations of the levels are found to differ 
(with respect to the standard computation without detailed 
follow-up of Z-states) by more than 20% for low redshifts 
[z < 800). In addition, the amplitude of the lines is also 
affected. Following all the angular momentum substates, 
the relative ratio of the different lines within a given se- 
ries increases. This can be seen, for example, in the case of 
the Balmer series, where the amplitude of the a-transition 
(An = 1) increases while the rest of transitions decrease. 



This shows that the simplifying assumption of statistical 
equilibrium within a given shell leads to an overall underes- 
timation of the net emission at high frequencies. 

We have also quantified the effects on the residual elec- 
tron fraction, showing that at low redshifts (z ~ 800) the 
x e is increased by 1.5%. The impact of this change on the 
angular power spectrum is small, producing differences of 
the order of 1% in the GYs. 

Finally, our computations have shown that, if the pri- 
mordial radiation field is assumed to be a pure black- 
body spectrum, then there is no significant emission from 
any line for redshifts greater than z = 2000, so a "pre- 
recombination" peak does not exists. 
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